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The extension of quantum trajectory theory to incorporate realistic imperfections in the measure- 
ment of solid-state qubits is important for quantum computation, particularly for the purposes of 
state preparation and error- correction as well as for readout of computations. Previously this has 
been achieved for low-frequency (dc) weak measurements. In this paper we extend realistic quan- 
tum trajectory theory to include radio frequency (rf ) weak measurements where a low-transparency 
quantum point contact (QPC), coupled to a charge qubit, is used to damp a classical oscillator cir- 
cuit. The resulting realistic quantum trajectory equation must be solved numerically. We present 
an analytical result for the limit of large dissipation within the oscillator (relative to the QPC), 
where the oscillator slaves to the qubit. The rf+dc mode of operation is considered. Here the QPC 
is biased (dc) as well as subjected to a small-amplitude sinusoidal carrier signal (rf). The rf+dc 
QPC is shown to be a low-efficiency charge qubit detector, that may nevertheless be higher than 
the dc-QPC (which is subject to 1// noise). 

PACS numbers: 73.23.Hk, 03.67.Lx 
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I. INTRODUCTION 



Solid-state proposals for building scalable quantum 
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seem promising. 



In 



information processor; 
any quantum information processor, the quantum bits 
(qubits) of information need to be read out, as well 
as controlled (via measurement-based feedback, for ex- 
ample). Quantum trajectory theorj^Z.^ has been used 
to describe single realizations of the continuous-in- 
time weak quantum measurement of electronic charge 
qubit^^^iH^^^ii^^ conditioned by the elec- 
trical output of a mesoscopic measurement device such as 
a quantum point contact (QPC) or single-electron tran- 
sistor (SET). Very recently, quantum trajectory theory 
has been applied to circuit QEDi^S. In these works the 
qubit evolution was conditioned on idealized measure- 
ment results (such as electron tunnelling) at the scale of 
the mesoscopic detector. That is, the observer's state 
of knowledge about the qubit state is updated based on 
these idealized measurement results. Two of us recently 
extended such work to condition the qubit state on a 
macroscopic signal that is realistically available to an ob- 
server in a dc-QPC measurement 3i In particular, this 
extension considered the noisy, filtering characteristic of 
an external circuit, including an amplifier. The result is 
a corrupted version of the idealized measurement results 
upon which the qubit evolution can be conditioned. This 
extension is known as "realistic quantum trajectory" the- 
ory and was pioneered for photodetection in quantum 
optics ; 22 ' 23 where it was applied to an avalanche photo- 



diode and a photoreceiver. 

In traditional dc charge-qubit measurement tech- 
niques, low frequency noise ("1// noise" 24 ' 25 ) limits the 
detector sensitivity^ To circumvent this, Schoclkopf 
et al^ introduced the so-called radio-frequency single- 
electron transistor (rf-SET). For the original configura- 
tion, the rf-SET demonstrated constant gain from dc to 
100MHz — an improvement on the conventional SET 
bandwidth by two orders of magnitude^ The idea is 
to measure the damping of a resonant (oscillator) cir- 
cuit in which the SET is embedded. In the context of 
charge-qubit detection^ 2 -! the damping depends on the 
qubit state, via the SET. Thus, monitoring the damping 
of the oscillator constitutes a continuous-in-time mea- 
surement of the charge qubit. This concept can be ap- 
plied to any charge-sensitive detector, in particular to the 
QPC, for example^ 

In this paper we derive an evolution equation for the 
conditional state of a charge qubit monitored continu- 
ously in time by a detector operated in the rf configu- 
ration. To the best of our knowledge, an equation of 
this type has not previously been derived. We also con- 
sider conditioning the qubit state on measurement results 
available to a realistic observer, within the framework 
of realistic quantum trajectory theory! 21 ' 22 ' 23 ' 29 In this 
approach, the bare charge-qubit detector (QPC in our 
case) is embedded in a realistic circuit, and an equation 
is derived that describes the evolution of the combined 
circuit-plus-qubit state conditioned on measurement re- 
sults available to a realistic observer. 
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Being able to determine the state of a quantum sys- 
tem conditioned on actual measurement results is ex- 
pected to be vitally important for quantum computa- 
tion, particularly for state preparation and quantum er- 
ror correction i 30 ' 31 i 32 ' 33 More broadly, it is also essen- 
tial for understanding and designing optimal feedback 

Control . !0-34. 35. 36. 37.38. 39.40. 41. 42. 43 

To simplify our analysis, we make a number of approx- 
imations. First, we use the low-transparency QPC as 
the charge-sensitive detector embedded in the rf circuit, 
instead of the SET. Second, we make a rotating wave 
approximation (RWA) to remove the 1// noise from our 
equations (as the rf configuration removes it in the exper- 
iment). In order to do this, we assume the QPC to be op- 
erating in the weakly- responding (diffusive) limiti^ ' 10 i 12 
In this limit, the QPC shot noise appears as white noise, 
which is modeled as a Gaussian-distributed stochastic 
process (having a diffusive appearance). Third, we as- 
sume that the rf-QPC is operated in the rf+dc mode 
introduced in the context of a SET in Ref. 0. This is 
where the QPC is subjected to a small amplitude sinu- 
soidal oscillation (rf ) superposed on a relatively large bias 
(dc). We find that the rf+dc QPC is a highly inefficient 
chargc-qubit detector. However, this may be higher than 
the measurement efficiency of the dc-QPC, which in prac- 
tice is not only limited by 1/ f noise, but is also further 
degraded the longer the measurement has to continue. 

The paper is organized as follows. The qubit and QPC 
are discussed in Sec. [TTJ The classical circuit is discussed 
and analyzed in Sec. IIII1 including the presentation of 
stochastic differential equations describing the state of 
the oscillator circuit. The stochastic master equation for 
the qubit state conditioned on the bare detector output 
is presented in Sec. [IVJ The realistic quantum trajec- 
tory derivation then proceeds in an analogous manner to 
Ref. Hi], where the circuit state was described by one pa- 
rameter — the charge on a capacitor. In general the oscil- 
lator circuit needs two parameters to describe its state, 
thus hinting at a complication in deriving the realistic 
equations for the rf configuration. However, because 
the QPC only damps the oscillator circuit (and doesn't 
induce a phase-shift in the reflected signal), all of the 
qubit information is contained in the amplitude (damp- 
ing) of the reflected signal. Therefore we can proceed 
with the realistic quantum trajectory derivation using 
a single parameter to describe the state of the classi- 
cal circuit. As always, numerical calculations must be 
performed to obtain the realistic quantum trajectories. 
However, in Sec. I VII we present analytical results for the 
heavily damped limit where the circuit is slaved to the 
qubit (adiabatic elimination of the circuit) . The paper is 
concluded in Sec. IVIII 



II. QUBIT AND QPC 

The measured quantum system we consider is 
the double-quantum-dot (DQD) charge qubit. 45 i 46 A 



schematic of the isolated DQD 4 ^ and capacitively cou- 
pled QPC is shown in Fig. [TJ We consider the low- 
transparency QPC, and represent it as a tunnel barrier 
between source and drain leads (reservoirs) with respec- 
tive Fermi levels /xs and //d- The QPC voltage bias is 
eVd = /ig — [ID, where e > is the quantum of electronic 
charged The DQDs are occupied by a single excess elec- 
tron, the location of which determines the charge state of 
the qubit. The charge basis states are denoted |0) and |1) 
(see Fig.[T]). We assume that each quantum dot has only 
one single-electron energy level available for occupation 
by the qubit electron, denoted by E\ and Eq for the near 
and far dot, respectively. 



DQD qubit 





Quantum point contact (QPC) 



FIG. 1: Schematic of an isolated DQD qubit and capacitively 
coupled low-transparency QPC between source (S) and drain 
(D) leads. 

The Hamiltonian for the qubit can be written as 



77, 



1 



qb 



(1) 



where the qubit energy asymmetry is e = E\ — Eq, Otun 
is the DQD tunnel-coupling strength, and a XlZ are Pauli 
matrices in the measurement (charge) basis. The eigen- 
values of H q h are +0/2, where SI = y/ Sl 2 un + e 2 . 

Associated with each of the qubit charge states is a cur- 
rent through the detector. The average current through 
the detector is I\ = eTi when the nearby dot is occu- 
pied, and Iq = eFo when the nearby dot is unoccupied. 
The variation in the detector output that depends on the 
qubit state is referred to as the detector's response and 
is denoted A7 = I\ — Iq. We can quantify the strength 
of the detector response by < |A7|/7 av < 2, where 
7 av = {1 1 + Iq) /2. Thus, a weakly responding^ detector 
satisfies |A7| <C 7 av , and a detector with finite, or strong, 
response satisfies |A7| ~ 7 av . In this paper we consider 
the limit of a weakly responding QPC, where the QPC 
shot noise appears as a diffusive, white noise process. 

The total Hamiltonian for the system is 
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Tot 



77 



qb 



i7p + 77r 



77, 



(2) 



where the qubit Hamiltonian 7f q b is given by Eq. (fT]). 
The free Hamiltonian describing the continua of electron 
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channels (momenta) k and q in the source and drain leads 



is 



H F = Y, ^sfc a L a sfc + w D g ar) g «D«. ( 3 ) 



where as and Sd are the Fermi field annihilation oper- 
ators for the source and drain leads, respectively. The 
tunnelling Hamiltonian 



Ht = ^2 T k q & Sk & Uq + T 9 *fe4,« 

k,q 



S/,- 



describes tunnelling between the source and drain leads. 
The probability amplitude for a source electron in chan- 
nel k to tunnel through the QPC into the drain channel 
q is T kq . The coupling Hamiltonian 



H, 



coup 



J2x kq & lk d 'D q +X* q k d 'u q d 'Sk 



(5) 



k,q 



describes the change in the effective QPC tunnelling am- 
plitude from Tkq — * Tkq + Xkq when the nearby dot is oc- 
cupied. This changes the QPC current from Iq = e |T| 2 
to ii = e\T + X\ 2 . Here T oc T^ q and X oc \kq are 
both proportional to the square-root of the source-drain 
voltage Vd-— ^ The occupation number operator of the 
nearby dot is h = (1 + a z )/2. Note that the height of 
the QPC (tunnel-junction) barrier is increased when the 
nearby dot is occupied, due to electrostatic repulsion, so 
that Iq > I\. 



III. CLASSICAL SYSTEM: OSCILLATOR 

An oscillator circuit, or tank circuit, consists of an in- 
ductor L, and capacitor C. We treat the oscillator clas- 
sically. The (angular) frequency for which resonance oc- 
curs in such an unloaded tank circuit, cjq = l/\/LC, is 
known as the resonance frequency. Embedding a dis- 
sipative component, like a resistor, into the tank cir- 
cuit provides a source of damping. For our purposes, 
the QPC provides the damping, so the oscillator damp- 
ing depends on the qubit state via the qubit-dependent 
QPC conductance. It therefore makes sense to moni- 
tor the damping of the tank circuit in order to ascertain 
the qubit state. This is achieved by using the tank cir- 
cuit to terminate a transmission line. Impedance mis- 
match between the tank circuit and transmission line 
causes a signal launched towards the tank circuit to be 
reflected back along the transmission line, where the re- 
flected signal can be observed^ See Fig.[3]for our equiv- 
alent circuit representing the rf-QPC setup. A recent 
experiment^ reports realization of a semiconductor rf- 
PC (radio-frequency point contact), with two benefits 
over the SET — lower PC impedance (thus simplify- 
ing impedance matching with the transmission line), and 



easier fabrication. This as yet unoptimized device ex- 
hibits a lower charge sensitivity than the rf-SET. In said 
experiment, the point contact is operating as a simple 
voltage-controlled resistor rather than a QPC. Another 
more recent experiment^ has realized fast charge sensing 
with a semiconductor rf-QPC. 

The AC voltage signal launched towards the tank cir- 
cuit is referred to as the "carrier" ^ Using a carrier sig- 
nal frequency equal to the resonance frequency of the 
unloaded tank circuit allows the detector to be replaced 

44,48 



(4) by its (frequency-dependent) dynamic resistance i?d 



This is a valid first-order approximation since the tank 
circuit is most sensitive to frequencies within a small 
bandwidth around the resonance frequency. 

Consider the equivalent circuit of Fig. [51 The oscillator 
circuit consisting of an inductance L and capacitance C 
terminates the transmission line of impedance Ztl = 
5051. The voltages (potential drops) across the oscillator 
components can be written as 



Vb(t) 
Vc(t) 



*(*), 
Qit) 
C ' 



(6a) 
(6b) 



where $(t) is the flux through the inductor, and Q(t) is 
the charge on the capacitor. These variables represent 
conjugate variables that together completely character- 
ize the classical oscillator state. The voltage across the 
detector, Vd(t), is the topic of the following subsection. 
The current flowing through the detector, I(t), will be 
discussed in Sec. [IVJ We define the incoming (relative to 
the tank circuit) transmission line voltage as the carrier 
signal V c (t) = Vjjj cos(wo^) plus some noise N(t): 



V in (t) = V c (t) + N(t). 



(7) 



The outgoing voltage in the transmission line is denoted 
Vout(t). As in the experiments of Ref. |49|, the detector 
bias Vq is added at a bias tee at the end of the transmis- 
sion line. 



A. Detector voltage 

In Sec. IIV1 we will discuss the conditional dynamics 
of the quantum system undergoing rf-QPC monitoring. 
Since the QPC tunnel ling rates depend on the detector 
voltage Vd(t) (see Ref. Qj]), we calculate it here. 

Using complex phasor notation, the voltage divider 
rule gives the detector voltage as 



V d (t) = V Q + Re 



Z 



V{t) 



since the inductor is a 



"short-circuit'' 

-l 



(8) 



at dc. The 



impedances here are Z^q = R^ L + iuj^C (i?dC combi- 
nation), and Z = Zb,c + (entire tank circuit), with 
Zl = iujqL (inductor). Here the AC voltage in the trans- 
mission line, 



V(t) = (l + (3)V c (t), 



(9) 
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FIG. 2: Equivalent circuit for continuous monitoring of a 
charge qubit coupled to a classical LC oscillator with induc- 
tance L and capacitance C. We consider the charge-sensitive 
detector that loads the oscillator circuit to be a QPC (see Fig. 
[l]for details). Measurement is achieved using reflection with 
the input voltage, Vi n (i), and the output voltage, V ut(t), be- 
ing separated by a directional coupler. The output voltage 
is then amplified and mixed with a local oscillator, LO, and 
then measured. 



consists of the incident and reflected AC signal. Here 
we make an assumption that N(t) comprises predom- 
inantly low-frequency noise which will be removed by 
the rf-QPC, and so we drop it to simplify the analy- 
sis. The reflection coefficient, (3, is given in terms of the 
impedance mismatch between tank circuit and transmis- 
sion line as 



Z + Ztl 



(10) 



That is, the AC signal reflected off the loaded tank circuit 
is given by (3V c (t). Combining these results, the detector 
voltage can be written as 



Vd(t) =^o + 2Re 



z 



RC 



Z + Z^l 
2Re [Ge**Vc(t)] . 



V c (t) 



(11) 



Here the amplitude gain G and phase-shift <f> of the carrier 
signal when it reaches the QPC are given in terms of the 
circuit and detector quality factors as 



G={Q~ d 2 Q- 2 



-1/2 



tan ((f)) 



(12a) 
(12b) 



The quality factors of the unloaded tank circuit, detector, 
and loaded tank circuit are Q = W0/7, Q& = 7d/wo, an< l 
Qt 1 = Q 1 + Qd > respectively. 

In the high quality limit Qd, Q > 1 (7d > w > 7), 
the tank circuit damping is due primarily to the detector. 
This limit therefore represents the highest sensitivity for 
the rf+dc QPC. We find that G ~ Qt and <j> ~ tt/2, so 
that the voltage across the QPC experiences a negative 
7r/2 phase-shift (cos — > sin) after the inductor. This is 
expected since the voltage across an (ideal) inductor lags 
the driving signal by 7r/2. In this limit the voltage across 
the QPC is 



V d (t) =Vo[l + £inSin(a;oi)]. 



(13) 



Here we have defined the dimcnsionlcss parameter e in = 
2QtV [ t J /Vq, which satisfies e- m <C 1 in the rf+dc mode4^ 
Equation (fTo^ is in agreement with Ref. |48| . 



B. Idealized classical dynamics 

The two conjugate parameters we use to describe the 
oscillator state are the flux through the inductor, $(£), 
and the charge on the capacitor, Q(t). The dynamics 
of the oscillator are found by analyzing the equivalent 
circuit of Fig. [5] using the well-known Kirchhoff circuit 
laws. Doing this we find that the classical system obeys 
the following set of coupled differential equations 



Q(t) 



-7$(t) 



Q(t) 
c 



V -2[V c (t) + N(t)}, 

(14a) 



V oat (t) = V c (t)+N(t) 



Z Th $(t) 



(14b) 
(14c) 



where I(i) is detector current and we have re-included the 
input noise N(t). These are the equations of a damped 
harmonic oscillator driven by both the input voltage and 
the QPC current (including the shot noise) at frequency 
ljq. The 1// noise in the input signal and the QPC is 
filtered out. To see this we start by recasting the problem 
in terms of the following dimcnsionlcss parameters 



it) 



1 



hZ LC 

Zlc 



[*(*) - $ ss ] = a x [*(*) - $ ss ] , (15a) 



Q(t) - Q ss ] = a y [Q(t) - Q ss ] , (15b) 



where a x and a y are implicitly defined above and the 
subscript ss refers to the steady state. The resulting 
equations are 

x(t) = -jx(t) + L0 y(t) - 2a x [V c (t) + N(t)}, (16a) 
y(t) = -w x(t) - a y I(t), (16b) 
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where I(t) = I(t) — I ss . The solution of these coupled 
equations involves the two time-scales uj^ 1 and 7 _1 . In 
the limit 7 _1 ^> ujq 1 (Q ^> 1), we can define a coarse- 
graining^ time-scale dt (Roman font d) that is short 
compared to 7 _1 , and long compared to uJq 1 . On this 
time-scale we can make the standard rotating wave ap- 
proximation (RWA), this being terms oscillating at fre- 
quencies greater than or equal to loq take their average 
value (unless they multiply white noise terms, that have 
non-negligible components at loq). 

Applying the above to Eqs. (flU)) gives 



dx(t) 

am 



ttx^in + %J 



(0 



dt 



a y I c ° s 



(t) 



dt, 



(17a) 
(17b) 



where the tilde denotes that we are in the frame rotating 
at ujq and we have dropped the effect of the input noise, 
N(t), as we assume that it is mainly 1/f noise and has 
negligible spectral weight at ui . The two currents I cos (t) 
and I sln (t) are the two quadratures of the filtered detec- 
tor current I(t) centered at the frequency loq. That is, 
they are given by 



T 



I s 



s dl 



I(s) cos(w s)ds, (18a) 



t-dt 
t 



I(s) sin(w s)ds. (18b) 



These components are band-pass-filtered currents over 
the bandwidth dt^ 1 and are the quantities to which an 
idealized observer would have access. To be more specific, 
the current coming from the QPC is 



I(t) =e\T\J(t)+I clas (t), 
where J(t) is the quantum signal given by 



j(t) = v«oW(**) +£(*)■ 



(19) 



(20) 



where the rate kq will be discussed later, and is the 
QPC shot noise, which we take to be white in the limit 
of a weakly responding QPCM^ That is, the QPC shot 
noise is Gaussian noise with the following correlations 



EK(t)] = 

mm(t')] = s(t-t'). 



(21) 



Here E denotes an ensemble average over all possible re- 
alizations of £(f). In Eq. (frnj) I c \as(t) is the deterministic 
classical component of the current (with the steady-state 
current subtracted — see Appendix [X] for details). Using 
Eq. (fl3|) we can write 



\T{t)\ « \T \ [1 + ±e in sin (w„t)] 



(22) 



since e ni <C 1 (see Sec lIIIA]) . Using this, the two quadra- 
ture currents become 



JCOS 

jsin 



(*) 
(t) 



e\T \VS^J cos (t) 
e|To|Vs^\/ sin (i)- 



tcos 
clas 



(*), 



(23a) 
(23b) 



where J cos (t) and J sln (t) are the two quadrature compo- 
nents of J(t). They are 



J COB {t) = 
J shl (t) = 



1 



i 



(24a) 



(24b) 



where S cos and S sln are dimensionless and defined such 
that [J cos (t)dt] 2 = [J sin (t)dt] 2 = dt. The Fourier compo- 
nents of the quantum signal J(t)dt are 



•C n (*) 



1 

dl 
i 

at j t _ 



J(s) cos(nw s)ds, (25a) 



t-dt 

t 



J(s) sin(nujQs)ds. (25b) 



-<i/ 



Note that [J^ (t)dt] 2 = [JZimt? = at/2 for n = 1,2; 
and that [J$ os (t)dt] 2 = dt [using Eq. (|2T|)]. 

To measure the two quadrature currents in Eq. (|23|) 
[and hence the quantum signals J cos (t) and J sin (t)], an 
ideal observer would measure the desired quadrature of 
the output voltage. This can be done by beating the 
output voltage with a local oscillator^ Using Eq. (|14c[) . 
the two output voltage quadratures are 



Kut(*) 

where B^„ 



^1 f ,r-( s ) e -^)/ 2 ds, (26a) 
Qx 2 J 



■Bcos 7 
«x 2 J 

e\T \a y 



jco S ( s)e -7(*-s)/2 dS; ( 26b ) 



and B r 



e\T \a y vS^ 



have units of t 1 / 2 . That is, B 2 in is the proportionality 
constant that turns the quantum signal, J(t), into the 
dimensionless current a y I mn , which drives the classical 
circuit (the same is also true for the cosine quadrature 
component). Here we have done some post-processing, 
removing the deterministic parts of V^(t) and V^(t), 
and dropping all contributions from N(t) as again we 
assume that it has negligible spectral weight at uj . By 
inverting these convolutions, the idealized observer can 
access both J cos (t) and J sln (i), enabling them to condi- 
tion the quantum state on these currents. The resulting 
equation describing the idealized conditional dynamics 
of the qubit is called a stochastic master equation, or 
quantum trajectory equation. Note that we can write 
Eq. ()26a|) as V™t(t) = jx(t) / (2a x ) (after some post pro- 
cessing), and similarly for Eq. (|26b[) in terms of y. 



IV. IDEALIZED QUANTUM DYNAMICS 

We now consider the idealized case where the stochas- 
tic QPC current can be perfectly measured [we have ac- 
cess to both J cos (t) and J sin (t)]. To describe the ideal- 
ized conditional qubit dynamics, we start with the mi- 
croscopic model of charge qubit monitoring by a (dc) 
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QPC in Ref. [l2|. We make an important modification 
to the model that is due to the time-dependence of the 
voltage across the rf-QPC, V d (t) [sec Eq. (jl3|) ]. It re- 
sults in time-dependent QPC tunnelling rates. We also 
make a rotating wave approximation (RWA) to simplify 
the analysis. The RWA is only possible for a weakly re- 
sponding QPC. In this paper, as in Ref. \l2\ we will refer 
to the limit of weak response as quantum diffusion. This 
is because in this limit there are many electrons passing 
through the QPC with each containing only a little in- 
formation about the qubit state, which under monitoring 
of the QPC current the evolution of the qubit will slowly 
wander towards one of the a z eigenstates rather then a 
sudden collapse. The linear form of the qubit stochastic 
master equation of Ref. [H is (see Appendix lAlfor details) 



dp c (t) 



H' b (t)+Hj(t),p c (t) 



dt 



|r d (t)D [a z }p c {t)dt 



+ [J(t)-fj]dtH ^{t)a z /2 p c (t), (27) 



where 



Kbit) = H qb + ha z \T(t)\ \X{t)\ sin(0)/2, 



Hj(t) = -M z [J{t) - n] vV(t)/2, 



(28a) 
(28b) 



and p is the mean of the ostensible distribution used to 
derive Eq. (|2"T)) (see Appendix [A]) . The linear supcropcr- 
ator, 7i, in Eq. (|2T[) is defined for arbitrary operators c 
by 

at 



H[c]p = cp + pc* 



pp. 



(29) 



It represents the back action effects of the continuous 
measurement. The timc-dcpcndcnt qubit dephasing rate 
is Td(i) = [^o(i) + Ki(t)]/2, where Ko(t) and Ki(t) repre- 
sent two types of measurement-induced dephasing in the 
qubit: Ko(t) represents information-limited (Heisenberg) 
dephasing^ which reflects the quantum-mechanical cost 
of obtaining information about the qubit state; K\{t) rep- 
resents additional (non-Hciscnbcrg) dephasing by pro- 
cesses that yield no qubit information. We define these 
dephasing rates by 



\Ao(*) = I x (t) I cos (0) = y/2T d (t) cos (9) , (30a) 

VM*) = \X(t)\sm(9) = y/2T d (t) sin (9) , (30b) 

where 9 is the relative phase between the QPC tunnelling 
amplitudes T and X. Note that if 9 = or 9 = tt, then 

K\ = 0, and Hj(t) = 0. This is a necessary condition for 
the QPC to be considered an ideal charge qubit detector 
operating at the quantum limit s 52 ' 53 
From Eq. (|13p we can write 



|*(t)l»v^[l 



rein sin 



(wo*)] , 



(31) 



since e m <C 1 (see Sec. lIII A)) . Here the time-dependence 
has been made explicit, and k/2 = \X\ 2 /2 is the time- 
independent dephasing rate for the charge qubit moni- 
tored by a dc-QPC (the situation in Ref. Il2l ). Substi- 
tuting Ko(t) and K\(t) (with the above approximation 



for their time dependence) into Eq. (|27[) and making the 
RWA gives 



dpc(t) 



d* 



H qh + Hj(t),p c (t) 
+ kV [a z ]p c {t)dt/A 
+ {j^(t)-p+ e f[JZ(t)-p]}dt 
x ft [^cos{9)a z /2] p c (t), (32) 

where J§° s {t) and J^(t) are defined in Eqs. (25]). The 
rotated versions of the Hamiltonians are 



1} 



H 3 (t) = ha z [Jr(t) - M + f K(t) A* 

x V«sin(0)/2, (33a) 
H qh ^H qh + ha z \T \^shi{9)/2. (33b) 



To get the linear form of the ideal quantum trajec- 
tory we need to rewrite Eq. (|32[) in terms of the signal 
that an ideal observer could access, namely J sm (i). This 
is achieved by expressing Jo° s (t) and J™(i) in terms of 
both the observed process, J sln (t), plu s some other un- 
observed processes, as done in Ref. [2l|. Averaging over 
the unobserved processes results in the following linear 
quantum trajectory equation in the RWA: 

dp c (t) = Cp c (t)dt 

+ VV [J sin (t)-Snfi] dtft p c (t), (34) 

where we have defined the efficiency r\ by 



V = 4 



t cos 2 (9)S S 
and the Liouvillian C is 



t cos 2 (0)/4, 



(35) 



Cp 



H qh + H'j(t) 7 p +-V[a z ]p. (36) 



We now have H' } (t) = - 
tion (l34l) normalizes to 



ha z J sin (t)yffjEtan{9)/2. Equa- 



d Pc (t) = C Pc dt + 07 [J s ' m (t) - y/rjK{&*)] dt 

x H [V^a z /2] p c (t), (37) 

where the normalized current is 



J sin (t)dt = y/TjR (& z ) dt + dW{t) 



(38) 



where dW(t) = £(t)dt is a Wiener increment^ This 
is the explicit expression for the quantum signal to 
which an ideal observer would have direct access. To 
get the correct statistics, we have made a Girsanov 
transformation,— which results in replacing p in the 
above with tJk(& z ). The stochastic Hamiltonian H' 3 {f) 
in Eq. (|36[) will also be updated in the same manner (re- 
placing p with y / k(it z )). The nonlinear back action su- 
peroperator is defined by its action on p with an arbitrary 
operator c by 



H[c]p = cp + pc) — (c + & ) p. 



(39) 
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Equation (|37|) is the first stochastic master equation 
presented for continuous measurement of a charge qubit 
using an rf configuration. It reveals some interesting 
physics about the rf-QPC operating in the rf+dc mode. 
First, we note that the qubit dephasing is due only to the 
dc component of the detector voltage, since k is a func- 
tion of Vq [and is independent of the AC component of 
Vd(i)]. Second, the rf-QPC is a highly inefficient detec- 
tor when operated in the rf+dc mode since rj oc e 2 n <C 1. 
Physically, this low detection efficiency arises because the 
qubit information is extracted relatively slowly by the 
(small) rf component of the QPC voltage, compared to 
the qubit dephasing by the (large) dc component. We 
conjecture that the low measurement efficiency of the 
rf+dc QPC may in practice be higher than that of the 
dc-QPC for two reasons: 1// noise limits the dc-QPC 
measurement efficiency and also further degrades the ef- 
ficiency the longer the measurement has to continue. 



V. REALISTIC DYNAMICS: REALISTIC 
QUANTUM TRAJECTORY EQUATION 

The previous sections conditioned the qubit evolution 
and the classical oscillator evolution on idealized mea- 
surement results available only to a hypothetical ob- 
server. For an experimentalist, it is much more useful to 
consider how to describe the qubit evolution conditioned 
on measurement results available in the laboratory. In 
rf-QPC or rf-SET experiments, the voltage leaving the 
transmission line is observed using homodyne detection 
of the amplitude quadrature {x ) . The phase quadrature 
(y) can be ignored because it is independent of the qubit- 
information-carrying signal, J sln (t) [sec Eqs. (fTT]) and 

mi 

In simple dyne detection (see schematic in Fig. [5]), the 
output signal V ou t(t) is amplified, and mixed with a local 
oscillator (LO). The LO for homodyne detection of the 
amplitude quadrature is Vlo (t) cos (u>ot) , where the 
LO frequency is the same as the signal of interest (or very 
slightly detuned). The resulting low-frequency beats due 
to mixing the signal with the LO are easily detected. 

The signal resulting from the homodyne detection is 



V{t)dt = y/A : 



SN 



■V c ™(x)dt + dW ont {t) 



\ s dt + dW out (t), 



(40) 



where Asn is the ratio of gain squared to noise power 
which has the dimensions of inverse time. This relates to 
the dimensionless signal-to-noise ratio for the measure- 
ment as SNR = As^t (t is the measurement duration). 
Equation (f4TH) shows that white noise is added to the 
amplified quadrature signal. Here 



SN 



7X 

2-Bsir 



(41) 



and the output noise Wiener increment dWo U t (t) satisfies 
[dWout(t)] 2 = dt. The addition of this extra noise pre- 



vents one from inverting the convolutions in Eqs. (|26[) 
to find J sln (i), and so the realistic quantum trajectory 
approac h 22 ' 23 must be employed to condition the qubit 
evolution on the realistic measurement record V(t). 

Our derivation of the realistic quantum trajectory 
equation describing the conditional evolution of the com- 
bined oscillator-qubit state closely follows the standard 
techniques ! 21 ' 23 The details are included in Appendix 151 
with the final result being the following superoperator 
Kushner-Stratonovich (SKS) equation: 



dp v (x) 



d_ 

dx 



BL^L 

2 dx 2 



dt 



+ (A2-A (2>p )[V(t)-A<2)Jdijpv(z) 
d 

- dt^fqRB sin — [a z p v (x) + p v {x)a z ] /2, 

(42) 

where we now consider x to be post-processed (the de- 
terministic part has been removed), and we have defined 
nix = —"fx/2. The true record the realistic observer 
would measure is 



V(t)dt = A <2) dt + dW{t), 



(43) 



where the observed white noise dW(t) is not the same 
as dW out (i) [contrast (03]) and (gOJ)]. In Eq. |g5J), the Li- 
ouvillian superoperator C contains the Hamiltonian evo- 
lution and qubit dephasing (and is defined in Appendix 
|B|) . The realistic quantum trajectory equation (|42|) is an 
important result of this paper. The first line describes 
the uncoupled, average evolution of the oscillator and 
the qubit. The second line describes the update of the 
realistic observer's state of knowledge of the supersystem, 
conditioned on the realistic homodyne output. The final 
line contains the effect of the qubit on the circuit. 

The (reduced) state of the qubit conditioned on the 
realistic output, py, is found from py(i) by integrating 
out the circuit variable x. Similarly, the (marginal) state 
of the circuit can be found by tracing out the qubit. 



VI. OSCILLATOR SLAVED TO QUBIT 

For sufficiently large damping 7 such that the first 
term in Eq. (|17ap dominates the dynamics of x (but still 
7d 3> 7, so that the RWA remains valid), the oscillator 
immediately damps to a qubit-dependent state. The cir- 
cuit is thus said to be slaved to the qubit. This adiabatic 
elimination of the circuit allows the conditioned qubit 
dynamics to be once again governed by the stochastic 
master equation alone. However, a realistic observer has 
access only to the homodyne output, that involves excess 
noise above that of J sln . Thus, we find the slaved output 
to which the realistic observer has access, then use it to 
condition the qubit state. The slaved value of x is found 
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by taking 7 — > 00 in Eq. (|26a[) . This allows us to make 
the standard delta-function replacement and write the 
quadrature output voltage as V^(t) = B s i n J sln (t)/a x . 
That is, the quadrature output voltage is directly pro- 
portional to the quantum signal an ideal observer would 
measure. Substituting this into Eq. ([40]) gives 

V sl (t)dt = ^/A^J sin (t)dt + dW out (t), (44) 

where J sln (i) is the signal used to condition the qubit 
state in Eq. ([34]) . Here we see that the output noise 
degrades the QPC signal, the effect of which is to reduce 
the efficiency of the detection. 

Using the above we can redefine the quantum sig- 
nal a realistic observer would measure as V s \(t) = 
J s i(i)V^SN + f , where 

J s i(t)dt = (&,) dt + dW(t), (45) 

where 



is the efficiency of the realistic rf+dc QPC in the slaved 
limit. Using this quantum signal, the stochastic master 
equation in the slaved limit is 

dp c {t) = dtCp c (t) + y/rjd [J s i(t) - y/rjsiK (a z )} dt 

x H [V^a z /2] Pc (t), (47) 

where the noise term in C is now in terms of J s i(t). The 
fact that the efficiency ?7 s i in Eq. (|4l)]) is strictly less than 
unity shows that the rf+dc QPC is incapable of reaching 
the quantum limit of a purity-preserving detector, even 
in the slaved limit. This conclusion is expected from 
the discussion following the idealized stochastic master 
equation in the RWA, Eq. ([57)1 . That is, even without 
including the circuit in our description of the rf-QPC, 
Eq. ([37]) shows that, for e- m <C f (the rf+dc mode), the 
rf-QPC is a highly inefficient charge qubit detector (that 
may nonetheless be more efficient than the dc-QPC due 
to the effects of 1/f noise in the dc case). Note that the 
theoretical maximum of r) s i is the idealized efficiency 77, 
which occurs for infinite signal-to-noise ratio in the am- 
plifier. As a final point, we note that the inclusion of 
input white noise in the analysis will have the effect of 
further reducing the rf+dc QPC efficiency (and signifi- 
cantly complicating the analysis)^ 



allows operation of the SET at high frequencies, where 
the 1/f noise prevalent in conventional low-frequency, dc 
measurements is completely negligible. 

We have two main results. We have derived an evolu- 
tion equation for the conditional state of a charge qubit 
monitored continuously in time by a detector operated 
in the rf configuration. To the best of our knowledge, an 
equation of this type has not previously been derived. 
This culminates in the stochastic master equation, or 
quantum trajectory, Eq. ([57)) . Our second main result is 
the extension of our idealized quantum trajectory (|37[) to 
consider conditioning the qubit state on a corrupted (fil- 
tered, more noisy) measurement signal available to a re- 
alistic observer. Our realistic quantum trajector y 21 i 22 i 23 
equation is Eq. (|42[) . 

Our model used the quantum-point contact (QPC) in 
place of the SET, and also assumed operation of the rf- 
QPC in the so-called^ rf+dc mode. In this mode, a small 
amplitude sinusoid (carrier) at the resonance frequency 
of the oscillator is superposed on a relatively large bias 
(dc). Experimentally^ this is done to maximize the sen- 
sitivity of a SET. In this mode we can make a rotating 
wave approximation (RWA) in the limit of a vanishingly 
small amplitude sinusoid. We found the rf+dc QPC to be 
a highly inefficient chargc-qubit detector, even ignoring 
the external circuit. Physically, this is because the small 
rf amplitude carries the qubit information, while the rel- 
atively large dc bias causes most of the measurement- 
induced qubit dephasing. It is important to note that 
the low measurement efficiency of the rf+dc QPC may 
in practice be higher than for the dc-QPC due to the 
debilitating effects of 1/f noise. Specifically, 1/f noise 
limits the dc-QPC measurement efficiency, and also fur- 
ther degrades the efficiency the longer the measurement 
has to continue. 

Having realistically modeled continuous-in-time qubit 
measurement using the rf+dc QPC, the next step is to 
extend realistic quantum trajectory theory to the pure rf 
mode (zero bias, larger rf amplitude) for both the QPC 
and the SET. These arc important future tasks, particu- 
larly for understanding and designing measurement and 
feedback control of quantum systems in the solid-state — 
valuable knowledge for fully harnessing the potential of 
future quantum technologies such as quantum comput- 
ers. 



VII. DISCUSSION AND SUMMARY 

In this paper we have presented a model for condi- 
tional monitoring of a charge qubit using an rf (radio- 
frequency) configuration. The rf configuration was intro- 
duced for the single-electron transistor (SET) in Ref. [26|. 
It involves embedding the SET in a resonant (oscilla- 
tor) circuit, and monitoring the resulting damping. This 
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APPENDIX A: RECASTING THE STOCHASTIC 
MASTER EQUATION OF GOAN, ET AL. 

The microscopic model of char ge q ubit monitoring by 
a conventional (dc) QPC in Ref. [t2| is used in this ap- 
pendix as a starting point for producing a stochastic 
master equation for the rf-QPC (in the rf+dc mode) [Eq. 
[27] . We add time dependence to the QPC tunnelling 
rates that is due to the time-dependent rf-QPC voltage 
(see Sec. lIII A]) , and recast the diffusive stochastic master 
equation of Ref. [l2| into a linear form. 

The diffusive stochastic master equation of Ref. UH is 



dp c (t) 



H q b,Pc(t) dt + V[T + Xn]p c (t)dt 



<{t)dt— [T*Xhp c {t) + X*T Pc (t)n] 
-t(t)dt^2Rc{T*X) (h) cPc (t) 



(Al) 



where n = |1)(1| and iJ q b, and the time-dependencies of 
T(t) and X(t) in our model are defined in Sec. [Ill This 
equation can be rewritten as 



dp c {t) 



Hqh(t),Pc{t) dt + |r d (t)D [&,] Pc (t)dt 



-n 



\/Ki(t)ia z /2 |p c (*), 



where 



H qh + ha z \T(t)\\X(t)\sm(6)/2, 

y/Kdfi=\X(t)\cOB(0), 

\X(t)\sm(e), 
/to(t)+Ki(t) \*(t)\ 2 



\/Kx{t) 

r d (*) 



(A2) 

(A3a) 
(A3b) 
(A3c) 

(A3d) 



Here 9 = axg(T*X) is the relative phase between T(t) 
and X{t). The superoperator H[c] in Eq. (|A2[) is defined 
by its action on p in Eq. ([39| . 

In Ref. |T3 the current through the QPC is 

I(t) = e{\T\ 2 + 2Re[T*X] (h) + \T\£(t)} , (A4) 

which in our parameterization is (now with time- 
dependent T and X) 



I(t)=e\T(t)\ [\T(t)\ + y/^t)(l + (a z ))+m\ -(A5) 

Here we see that this current comprises two parts: a 
quantum signal, J(t), that depends on the state of the 
qubit and the noise and a large deterministic clas- 
sical signal, I c \as(t)- The quantum and classical signals 
are 



J(t) = y/^tfj (a z ) + C(t), 



/clasW = e|T(t)| y/KoV) + \T(t)\ 



(A6) 
(A7) 



which allows us to write I(t) as 

I(t) = e\T(t)\J(t)+I clas (t). 



(A8) 



Using the quantum signal we can rewrite the stochastic 
master equation [Eq. (|A2[) ] as 

dpc(t) = -~ [fl^ b (t) + ffj(t) )/ 9 c (t)] d£ 

+ |r d (t)2?[CT«]p c (t)<ft 

+[J(t) - ^^t)(& z )}dtn \^K (t)a z /2^p c (t), 

(A9) 

where 



Hj(t) = -ha z J(t)- y/Ko® (a,) (A10) 



Equation (|A9[) is the normalized quantum trajectory an 
observer would use to describe their state of knowledge of 
the qubit if they had access to the quantum signal J{t). 

To derive the realistic quantum trajectory, we use lin- 
ear quantum trajectory theory^ i 29 i 57 i 58 so we now recast 
the above in linear form. The linear quantum trajectory 
is derived in essentially the same way as Eq. (|A1|) was 
done in Ref. [H, except that one must use an ostensible 
distribution for J rather then the real signal [Eq. IA6j . 
That is, the possible results Jdt at time t are chosen 
from a Gaussian distribution with mean p and variance 
dt. Here p is an arbitrary ostensible parameter, with the 
only constraint that the ostensible distribution for J must 
be non-zero when the true distribution is non-zero. We 
don't go through the calculation in detai l, but the final 
result is equivalent to simply replacing yf Ko(t) (a z ) with 
p in Eq. (|A9[) . That is, the linear quantum trajectory is 



dp c (t) 



dt 



H 3 (t),p c (t) 
f- \T A {t)V [&,] p c {t)dt 

-\j{t) - p}dtn [Vko®*,/ 2 ! 



(AH) 



where the linear measurement superoperator is 

H[c]p = cp+ p$ - pp, (A12) 

and 



H 3 (t) = -ha z [J(t) - fi] (A13) 
Equation ([ATT]) is Eq. ([27]) . 



APPENDIX B: DERIVATION OF THE 
REALISTIC QUANTUM TRAJECTORY 
EQUATION 

In this appendix we derive a realistic quantum tra- 
jectory equation for the rf-QPC-monitorcd charge qubit. 
The derivation closely follows the derivation presented in 



10 



Ref . l2ll for the conventional QPC. As discussed, we ignore 
the phase quadrature of the classical oscillator since it re- 
veals no qubit information [see Equations (|17[)]. There- 
fore, we describe the imperfect knowledge of the oscillator 
by a probability distribution for the amplitude quadra- 
ture, P(x). 



Stochastic Fokker-Planck equation for the 
oscillator 



For notational simplicity, we express the Langevin 
equation for x(t), Eq. (|17a[) (after the deterministic clas- 
sical part has been removed), as 



dx(t) = [ms + B sin J sin 



(*)] dt, 



where we have defined 



(Bl) 



(B2) 



Equation (|B1|) describes the evolution of x(t) for per- 
fect knowledge of J sul (t). A realistic observer will not 
have direct access to the idealized quadrature current in 
Eq. (IB II) . so we find an equation for P(x) (see Refs. [23l 
and l2ll for details of the procedure). The result is the 
stochastic Fokker-Planck equation: 

dP c (x) = dt j- A [ m£ + B siu J sin (t)] 



2 dx 2 



P(x), 



(B3) 



As expected, we have both deterministic and stochastic 
drift (d/dx) of P(x), as well as the diffusion {d 2 /dx 2 ) as- 
sociated with the stochastic drift. We next consider con- 
ditioning P(x) on the homodync measurement results. 



equation is arbitrary (subject to the condition that A(V) 
is nonzero when P(V) is nonzero). Inspection of Eq. (|4"0|) 
shows that P X {V) is a Gaussian distribution of mean Aj 
and variance v = l/dt. That is, 



dt 



/2n 



cxp ■ 



\v-\iYdt 



(B6) 



where Aj is defined in Eq. (|4Tj) . The results above com- 
bine into Eq. (|B4[) to give the Zakai equation 



(B7) 



P V {Z) = {1 + Vdt\ £ }P(x). 



This is exactly analogous to the Zakai equations derived 
for the conventional (dc-mode) QPC in Ref. l2ll . 



3. Combined equation for the oscillator 

As in Ref. I2H we choose to condition on measurement 
results after microscopic processes, so that the combined 
conditional evolution of the oscillator is given by 

P(x) + dP v , c (i) = {1 + VdtXi} [P(x) + dP c (5)] , (B8) 

where dP c (x) is given by the stochastic Fokker-Planck 
equation (|B3[) . We find that 



BL d 2 



dP v . c (x) = dt j - J, [ms + B sin J sin (t)] 



2 dx 2 

VA 5 \P{x). (B9) 



This equation describes the evolution of the oscillator 
state (via the amplitude quadrature) conditioned on both 
observed and unobserved processes. 



2. Zakai equation for the oscillator 

Following Ref. [U, we find the best estimate for 
P{x) conditioned upon the measurement result V using 
Baycsian analysis. Denoted P\;(x) 1 this estimate is 



Pv{Z) 



Px(V)P{x) 
A(V) : 



(B4) 



where the bar indicates an unnormalized distribution. 
The ostensible distribution A(V) is a Gaussian distribu- 
tion of arbitrary mean A, and some variance v: 



A(V) 



1 



: CXp 



(V-Af 
2v 



(B5) 



That is, for the Zakai equation we consider the observed 
output V to be ostensibly Gaussian white noise of mean 
A. We simplify the derivation by choosing A = 0, but 
note that the choice of ostensible mean (in A) in the Zakai 



4. Stochastic equation for the joint qubit-oscillator 

state 

The joint system (supersystem) evolution conditioned 
by all processes is found by joining the stochastic master 
equation [Eq. (|37p] and the combined classical evolution 
[Eq. (|B9[) ] . This step is performed vi a 21 ' 23 

Pv,c{x) + dp V:C {x) = [P ViC (x) + dPv, c (£)] 

x [p c (t) + dpc(f)] ■ (BIO) 

The coupled noise process in dPy jC [Eq. (|B9|) ] and dp c 
[Eq. I|34p] is J sln (t). This is the term from which corre- 
lations between the circuit and qubit evolution arise in 
our joint stochastic equation. This equation describes the 
evolution of the joint quantum-classical state conditioned 
on the observed process V(t) and the unobserved process 
J sln (i). To remove the conditioning on the unobserved 
process, we simply average over it. In our approach, the 
unobserved process is ostensibly a white noise of zero 
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mean (since we use the linear equation). The result is 
the following superoperator Zakai equation 



d 

dp v (x) = dtl - -Q^mx 



- dt,ffjKB sin — [a z p v (x) + p v (x)a z ] /2, 

(Bll) 

where 77 is defined in Eq. (|35p . and C is defined in 
Eq. ©. 

This description of the joint quantum-classical system 
evolution conditioned on the observed output V(t) does 
not preserve normalization of the state. This is because 
ostensible statistics were chosen for V(t). To complete 
the realistic quantum trajectory derivation, we must next 
find the true statistics for V(i), and normalize the super- 
operator Zakai equation. 



5. Superoperator Kushner-Stratonovich equation 

Normalization of the superoperator Zakai equation is 
performed^ by taking the trace over the qubit and inte- 
grating over the oscillator (a;): 



p v (x) +dp v (x) 



p v (x) + dp v {x) 
J Tr [pv(x) + dpy(i)] dx 



(B12) 



It is important to note that the resulting equation is not 
the superoperator Kushner-Stratonovich (SKS) equation 
because ostensible statistics are still being used for the 
realistic homodyne output V(t). To find the SKS equa- 
tion, we must substitute the true expression for Vdt. The 
true distribution for V is found from the superoperator 
Zakai equation (jBllj) by tracing over the qubit and in- 
tegrating over all x, then multiplying the result by the 
ostensible distribution A(V) [Eq. (|B5)) ]. That is, 



P(V) = A(V) / Tr [p v (x) + dpy(x)] dx 
Vdt 



exp 



V-A 



dt/2 



(B13) 



where A( 2 )^ = \Msn7 (x) p /2B s - ln [c.f. Eq. idU], and the 
average (x) p = J JTr [p{x)] dx is qubit-dependent (hence 
the p subscript) . This means that the true expression for 
V(t) is 



V(t)dt = A<2) dt + dW(t), 



(B14) 



where dW{t) is the white noise Wiener increment a re- 
alistic observer would see. The end result is the SKS 
equation 
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